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We suggest a technique to induce effective, controllable in- 
teractions between atoms that is based on Raman scattering 
into an optical mode propagating with a slow group veloc- 
ity. The resulting excitation corresponds to the creation of 
spin-flipped atomic pairs in a way that is analogous to cor- 
related photon emission in optical parametric amplification. 
The technique can be used for fast generation of entangled 
atomic ensembles, spin squeezing and applications in quan- 
tum information processing. 

PACS numbers 03.67.-a, 42.50.-p, 42.50.Gy 



The intriguing possibility for controlled manipulation 
of interacting quantum systems is the basis for a number 
of exciting developments in the field of quantum infor- 
mation science [y. These are expected to have an im- 
pact in a broad area ranging from quantum computation 
and quantum communication [Ej to precision measure- 
ments [Bj and controlled modeling of complex quantum 
phenomena |j]. 

This Letter describes a new technique to induce effec- 
tive coherent interactions between atoms in metastable 
states. The technique is based on a resonantly en- 
hanced nonlinear process involving Raman scattering 
into a "slow" optical mode p|, which creates a pair of 
spin-flipped atom and slowly propagating coupled exci- 
tation of light and matter (dark-state polariton) . When 
the group velocity of the polariton is reduced to zero [^,0 , 
this results in pairs of spin flipped atoms. 

The present phenomenon of spin pair creation exhibits 
strong similarities with optical parametric amplification 
(OPA), in which pairs of photons are generated that pos- 
sess non-classical correlations in photon number, quadra- 
ture component fluctuations or in polarization states [g| . 
In direct analogy, the present technique is capable of gen- 
erating non-classically correlated atomic ensembles and 
entangled spin excitations. The latter can easily be con- 
verted into corresponding states of photon wavepackets 
"on demand" , which makes the present approach most 
suitable for implementing protocols in quantum informa- 
tion processing that require a combination of determin- 
istic sources of entangled states and long-lived quantum 
memory ||9|,pl|. 

The present technique can also be viewed as a new 
mechanism for coherent "collisions" [11| between atoms 
mediated by light. In particular, the case when atomic 



pairs are excited into two different levels (as e.g. in 
Fig. la) closely resembles coherent spin-changing inter- 
actions that occur in degenerate atomic samples [E2I, 
whereas the case when atomic pairs are stimulated into 
identical states (Fig. lb) is reminiscent of dissociation of 
a molecular condensate Q . To put this analogy in per- 
spective we note that the rate of the present optically 
induced process can exceed that of weak interatomic in- 
teractions by orders of magnitude. Therefore the present 
work may open up interesting new possibilities for study- 
ing many-body phenomena of strongly interacting atoms. 

Before proceeding we note that a number of proposals 
have been made for generating entangled states of atomic 
ensembles and resulting in so-called spin squeezed states. 
Some are based on interatomic interactions at ultra-cold 
temperatures [[l4[, whereas others involve mapping the 
states of non-classical light fields into atoms [|5[, QND 
measurements of spins ||l^] with light and Rydberg block- 
ade |17|. In contrast to these mechanisms the present ap- 
proach does not require coherence of the atomic motion 
or sources of non-classical light and is completely deter- 
ministic thereby significantly simplifying possible exper- 
imental realizations. We further show that the present 
technique can be made robust with respect to realistic 
decoherence processes such as spontaneous emission and 
leakage of slow photons from the medium. We also note 
the work of Franson and co-workers pq | on quantum logic 
based on so-called "photon-exchange" interaction. How- 
ever further analysis fl^ showed that these mechanisms 
do not result in a non-linear interaction. 

We consider a system of N atoms (Fig.l) interacting 
with two classical driving fields and one quantized mode 
that is initially in a vacuum state. Relevant atomic sub- 
levels include two manifolds of metastable states (e.g hy- 
perfine sublevels of electronic ground state) and excited 
states that might be accessed by optical transitions. The 
atoms are initially prepared in their ground states \g). 
One of the classical fields (Rabi frequency fii) is detuned 
from the atomic resonance by an amount roughly equal 
to the frequency splitting between ground state mani- 
folds. The other (Rabi frequency ^2) is resonant with an 
atomic transition I62) -^ |a2). The quantized field can 
be involved in two Raman transitions corresponding to 
Stokes and anti-Stokes processes. Whereas the former 
corresponds to the usual Stokes scattering in the forward 
direction, the latter establishes an Electromagnetically 
Induced Transparency (EIT) and slows down its group 



velocity. The pair excitation can be viewed as resulting 
from quantized photon exchange between atoms (Fig. 2) 
in a two-step process. The first flipped spin is created 
due to Stokes Raman scattering, which also results in 
photon emission in a corresponding Stokes mode. In the 
presence of EIT, this photon is directly converted into a 
dark-state polariton which becomes purely atomic when 
the group velocity is reduced to zero. This implies that 
atomic spins are always flipped in pairs. In Fig. la the 
two final states involved in Raman transitions are differ- 
ent and atomic pairs in different states are created. In 
Fig. lb the final states of the two Raman processes are 
identical, in which case atomic pairs in the same state 
result. 

In what follows we will focus on a system (Fig. la) in- 
volving two atomic modes. Consideration of the scheme 
of Fig. lb proceeds along the same lines. For conceptual 
simplicity we here assume that the quantized field corre- 
sponds to a single mode of a running-wave cavity with 
a creation operator a^ and atom-field coupling constants 
gi and g2- Generalization to multi-mode i.e. travelling 
wave configuration is straightforward. The interaction 
Hamiltonian for the system of N atoms and light can be 
split into two parts H — Hram + Hres , which are given 
by: 
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[hg2a^^ga2 + hn2^b2a2 + h-c], 
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where T,^^ — J2^\^Aii{'^\ ^'^^ collective atomic opera- 
tors corresponding to transitions between atomic states 
\lj),\v), A is the detuning of the classical field fii from 
the single-photon transition \g) — > |ai), Si and 62 arc the 
two-photon dctunings from the \g) —f \bi) and \g) -^ I62) 
transitions respectively as shown in Fig.l. 

In the limit of large detuning A and ignoring two- 
photon detunings for the moment, the Hamiltonian 
Hram dcscribcs an off-resonant Raman scattering. Af- 
ter a canonical transformation corresponding to adiabatic 
elimination of the excited state Hram becomes equivalent 
to: 
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= hy/N^^aSi+ h.c. 
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where we disregarded the light shift Sl — il^/A and in- 
troduced ^i = l/\/]VSgbj . The light shift 6l can be eas- 
ily compensated by re-defining the energy of the atomic 
levels and will be disregarded in the remainder of this 
Letter. The resonant part of the Hamiltonian Hres is 
best analyzed in terms of dark and bright-state polari- 
tons 
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which are superpositions of photonic and atomic excita- 
tions, with S2 = l/vNT,gb2- In particular, Hres has an 
important family of dark-states: 
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with zero eigenenergies. Note that all other eigenstates 
of Hres have, in general, non-vanishing interaction en- 
ergy. Under conditions of Raman resonance and suffi- 
ciently slow excitation ( "adiabatic condition" ) the Stokes 
photons emitted by Raman scattering, Eq. (0) , will there- 
fore couple solely to the dark-states (|g). In this case the 
evolution of the entire system is described by an effective 
Hamiltonian: 



Heff = K{PdSi + sIpI), 



(6) 



with ^ = f^if^z/A X glVN/^/\g2\^N + nl. The Hamil- 
tonian (^) describes the coherent process of generation of 
pairs involving polaritons and spin-flipped atoms. Note 
that for a small number of excitations the spin waves and 
polaritons obey bosonic commutation relations and this 
Hamiltonian is formally equivalent to that describing op- 
tical parametric amplification (OPA) of two modes [||. 
This is also analogous to the "counter-twisting" model 
of Ref. ||21|], which is known to result in maximal spin 
squeezing for large number of excitations. 

We now consider the scenario in which the system is 
evolving for a time r under the Hamiltonian Heff , after 
which both fields are turned off. If the procedure is adia- 
batic upon turn-off of the coupling fields rii_2 the polari- 
tons are converted into pure spin excitations Pd -^ 82- 
Hence the entire procedure will correspond to the follow- 
ing state of the system: 
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Here |n(,j,nh2) = l/n!(S'^)"(S'J)"|g) are Dicke-like sym- 
metric states of atomic ensemble and we assumed 
nbi.b2 ^ N- For non-zero S,t this state describes an en- 
tangled state, for which relative fluctuations between the 
two modes decreases exponentially to values well below 
the standard quantum limit (SQL) corresponding to un- 
correlated atoms. 

The above analysis only includes the interaction with 
a single (forward-propagating) quantized radiation mode 
and neglects decoherence processes. We now take into 
account realistic decoherence mechanisms such as spon- 
taneous emission from the excited states in all directions 
and decay of the cavity mode with a rate k. The evolu- 
tion of atomic operators is then described by Heisenberg- 
Langevin equations: 



-'fiiy — 'yfiiy^fii^ \ , [-H , lu^y\ -r -T^i/, 
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where 7^1/ is a decay rate of coherence /i — > i^ and F^^, are 
associated noise forces. The latter have zero average and 
are ^-correlated with associated diffusion coefficients that 
can be found using the Einstein relations. The cavity 
mode obeys the equation of motion: 



d = -Kfl - igi'Eb^ai - i92^ga2 + Fa{t). 



(9) 



We proceed by adiabatic elimination of optical polar- 
izations associated with Stokes emission. To this end we 
assume large single-photon detuning A 3> 7 and to first 
order in a we obtain the following equations of motion 
for the metastable coherences 



(10) 



<5i = - [igb + i5i\ S[ + I — a + F'g^ 

82^- [igb + i52] S2 - l{n2/ ^)^ga, + Fg, 



where 7gb = ^gb+^L and 7^ = 7ag|Jlip/A2 is the optical 
pumping rate. The ground state coherence thus decays 
at a rate modified by isotropic spontaneous emission from 
the excited state (we consider the case 7l » Igb)- 

To treat the resonant EIT-like interaction we first 
rewrite the equations of motion in terms of the dark 
and bright polariton operators (^j and proceed to adi- 
abatically eliminate the optical coherence Sga2 and the 
bright state polariton Pg. In the relevant limit when 
\g2?N/^agH > 1 and when rj = |.g2piV/|Sl2p the ratio of 
vacuum light velocity to group velocity is large 77 ^ 1, 
we find: 

Pd = -\ ;— \Pd +^iS[ +FD[t), 



1+7? 



Sl = [g^7L -IL- t5i]Sl - liPn + Fl (t). (11) 

We note that cavity losses are strongly suppressed in 
the limit r? 3> 1: subsequent to the large group velocity 
reduction |q|, the polariton is almost purely atomic and 
the excitation leaks very slowly out of the medium. 

The equation of motion for coherence Sl contains a 
loss term (due to isotropic spontaneous emission) and 
a linear gain term (due to emission into bright polari- 
ton). The two can compensate each other. However the 
linear phase-insensitive amplification is also accompanied 
by correspondingly increased fluctuations, represented by 
new Langevin forces Fu (t) , Fg (t) . Terms resulting from 
the direct coupling of dark and bright polaritons can be 
neglected for |J72| S> JagK., which is reminiscent of the 
condition under which KIT is established. 

To quantify the resulting quantum correlations we 
introduce a squeezing parameter in direct analogy to 
the optical parametric case. We define the quadratures 
Xi = {Sl + Sl)/V2, Yi = i{Si - Sl)/V2; these can be 
measured e.g. by converting spin excitations to light. 



Correlations between the modes appear due to dynami- 
cal evolution and squeezing is found in the quadratures of 
the sum and difference modes X± = {Xi ± X-2)l \/2 and 
Y± — (Yi ± Y'2)\f^. For small number of excitations the 
sum and difference modes obey standard commutation 
relations [Xa,y^] = —i5a,p where a,(3 — +, — orl,2. A 
quadrature Y± is squeezed when Ay±(i)^ < 1/2. 

We find that squeezing is optimal under conditions of 
four-photon resonance {5i =82) and in the limit of 77 ^ 1 
(Fig. 3). Evolution leads to maximum squeezing of y+ at 
t = t* after which the growing fiuctuations in X+ give 
rise to increased noise in Y+ . Note that the number of ex- 
citations grows exponentially with time (Fig. 3c). Specif- 
ically, in the case gi — g2, for ^t > 1, we have: 

Ar+(t)2 = 1/2 |e-2e* + 2^ + ^ (12) 
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where we have neglected terms of higher order in jl/S. 
and k/^. The maximum amount of squeezing is ob- 
tained after an interaction time t* such that e~^^* = 
{iL + k/7?)/4^ and is given by AYf = {5jl + 3k/?7)/4,^. 
Since both the interaction parameter ^ and the relax- 
ation rate of the polariton "fjj =71,-1- k/tj depend on 
the single photon detuning A (Fig. 3a), we find that 
squeezing is optimized for the single-photon detuning 

Aopt - 7agA/ 3jn;p^;^;7-, and 



Ay; 



yis/i 



opt 



Vlg^FWi^ 
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Note that the denominator is equal to the atomic 
density-length product multiplied by an empty cavity 
finesse and can easily exceed 10^ even for modest val- 
ues of the density-length product and cavity finesse. 
Note also that the strong coupling regime of cavity QED 
g > max[K, 7] is not required to achieve strong correla- 
tions, in fact as long as g^N > 2kj squeezing is achieved. 
Furthermore, although a cavity configuration was used 
for simplicity, the results of the present analysis remain 
qualitatively valid in the limit of unity finesse, i.e. free 
space. We further emphasize that typical generation rate 
resulting in such optimal squeezing flifl2/A.opt can eas- 
ily be on the order of fraction of MHz. In such a case 
other decoherence mechanisms are negligible. Doppler 
shifts can also be disregarded as long as all fields are 
co-propagating. 

To summarize, we have presented a scheme based on 
the interaction of coherent classical light with an op- 
tically dense ensemble of atoms that leads to effective 
coherent spin-changing interactions involving pairs of 
atoms. We have shown that this process is robust with 
respect to realistic decoherence mechanisms and can re- 
sult in rapid generation of correlated (spin squeezed) 



atomic ensembles. Furthermore, the resulting spin ex- 
citations can be easily converted into photons on de- 
mand, which facilitates applications in quantum infor- 
mation processing. Possible applications involving high- 
precision measurements in atomic clocks can be also fore- 
seen. We further note that extension of this work into 
the domain of very large atomic density-length product 
or high-finesse cavities might allow to create maximally 
spin-squeezed states or macroscopic quantum superposi- 
tions ( "Schrodinger cat" states). This in turn might al- 
low to observe interaction-induced quantum phase tran- 
sitions [ p2[ |. 
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FIG. 1. Level scheme for the coherent interaction leading 
to pairs of atoms in (a) different final states [62) and |6i), (b) 
the same final state \b). 
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FIG. 2. Diagram illustrating coherent atom-atom interac- 
tion mediated by dark-state polariton, leading to the creation 
of a pair of spin-fiipped atoms. 
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FIG. 3. (a) Quadrature variance AY^ vs. single-photon 
detuning A and interaction time ^t, (b) same for A = Aopt 
and 5i = 82 showing maximum squeezing AY^. ~ 0.02 (for 
•\/t;|A/7agK = 100), (c) Number of excitations pumped in the 
system vs. time (same conditions as in b) and (d) AY+(i*)^ 
vs. two-photon detuning 5 = (5i — 52)/2 for A = Aopt and 
where i* gives maximum squeezing. 



